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Interrelated inequalities involving doubly stochastic matrices are presented. For example, if B 
is an n by n doubly stochastic matrix, x any nonnegative vector and y = Bx, then X\Xf • • •%% < 
J\j2' ' -Jn. Also, if A is an n by n nonnegative matrix and D and E are positive diagonal matrices 
such that B r: DAE is doubly stochastic, then det DE > p(A)~'\ where p(A) is the Perron- 
Frobenius eigenvalue of A. The relationship between these two inequalities is exhibited. 
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An n by n entry-wise nonnegative matrix B = (6,,) is called row (column) stochastic if 2 b\j = 1 

for all i = 1,' • ",ti ( 2 bij = 1 for all / = !,• • * 9 n J. If B is simultaneously row and column stochastic 

then B is said to be doubly stochastic. We shall denote the Perron-Frobenius (maximal) eigenvalue of 
an arbitrary n by n entry-wise nonnegative matrix A by p(A). Of course, if A is stochastic, p(A) = 1. 
It is known precisely which n by n nonnegative matrices may be diagonally scaled by positive 
diagonal matrices D, E so that 
(1) B = DAE 

is doubly stochastic. If there is such a pair /), E, we shall say that A has property (*). In this event it is 
our interest to obtain inequalities on D and E. In the process, certain related inequalities for doubly 
stochastic matrices are noticed. 

It was first realized by Sinkhorn [4] 1 that if A is entry-wise positive and square, then A has 
property ( * ) . The proof amounts to showing that the process of alternately scaling A to produce a row 
stochastic matrix, and then a column stochastic matrix, and then continuing the process, actually con- 
verges to a doubly stochastic matrix. The hypothesis of positivity, however, can be weakened somewhat. 
If there exists no single permutation matrix P such that 



pTAp =(t L) 



where An and A 2 2 are square, then A is called irreducible. If there is a pair of permutation matrices 
P, Q such that C = PAQ, then we shall say that A and C are equivalent. If, further, A is equivalent to 
no matrix of the form 
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(2) 



(An \ 



where An and A 22 are square, then A is termed completely irreducible. It is an easy calculation to show 
that A has property (*) if and only if each matrix equivalent to A does, and it is equally clear that for A 
to have property (*) it must have the zero-nonzero sign pattern of a doubly stochastic matrix. For 
example, if A has property (*) and if A is of form (2) this means we must then have A 2 i = 0. It has 
further been shown [1] that when A is completely irreducible the alternate scaling process of Sinkhorn 
still converges and thus A has property ( * ) . Since property ( * ) is preserved under direct summation, we 
may summarize as follows. 

Remark 1: A square nonnegative matrix A has property (*) if and only if A is completely irreducible 
or A is equivalent to a direct sum of completely irreducible matrices. 

Thus, property (*) depends only on the zero pattern of A. It is also a straightforward calculation (fol- 
lowing [4] ) that 

Remark 2; If A has property (*) , then the product DE of (1) is unique. 

Our first observation is both necessary for later proofs and of interest by itself. 

Theorem 1: // B = (by) is an n by n doubly stochastic matrix and x > is any nonnegative vector, 
then, for y = Bx, we have 

n n 

(3) IT x, < n y,. 

i=i i=i 

// B is completely irreducible, equality holds in (3) if and only if the right-hand side is or all com- 
ponents of x are the same. Furthermore, among all irreducible nonnegative square matrices B satisfying 
p(B) < 1, only those diagonally similar to doubly stochastic matrices satisfy (3) for all x > 0. 
Proof: From the arithmetic-geometric mean inequality [2] 



(4) 



n*i yi < 2 yiXi 



where x = (#i, # * ',x n ) T is any nonnegative vector, and y = (yi,* * • ,y n ) is a vector of nonnegative num- 

n 

bers satisfying 2 yt =■ 1. Equality holds in (4) if and only if the xfs corresponding to nonzero y/s are 
all equal. Now, suppose B = (fcy) is row stochastic and y = Bx, x > 0. It follows from (4) that 

n bij n 

(5) n Xj f < 2 bijXj = y^ for i = 1 5 * • *,n. 

;=1 ' " /=! 

Taking a product over i of both sides, we arrive at 

n n n 

(6) IT Xj2>bij < n y^ 

;=1 i=l i=\ 

n 

If B is doubly stochastic, 2 b^ = 1 for each / = 1, • • • ,n, and it follows that (3) holds. 

i=i 

To analyze the case of equality, it is clear that equality holds in (3) if either x is a vector of equal 
components or the right-hand side of (3) is 0. On the other hand, if equality holds in (3) and the right- 
hand side of (3) is not 0, then equality must hold in (5) for each i = 1,* • - ,n. This means that for each 
i, the #/s corresponding to nonzero 6j/s are all equal. This, in turn, implies, by virtue of equality hold- 
ing in (5) for all i, that y = Q T x for some permutation matrix Q. Since BQQ T x = y, we have that BQ 
has Q T x = y as a Perron-Frobenius eigenvector (corresponding to p(BQ) = I). If B is completely 
irreducible, then BQ is irreducible, and, since BQ is doubly stochastic, its Perron-Frobenius eigenspace 
is one-dimensional (because of the irreducibility) and is spanned by (1,1, • * *,1) T . Therefore all com- 
ponents of Q T x, and thus of x, are the same. It should be noted that even in case B is completely irreduc- 
ible it is possible that the right-hand side of (3) be for a nonnegative nonzero vector x. For example, let 
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1/3 1/3 1/3 \ / 

B = [ 1/6 1/6 2/3 and* = 

1/2 1/2 0/ \ 1 

Thus, the second statement of the theorem is proven. It should also be noted that the case of equality in 
(3) may similarly be analyzed if B is equivalent to a direct sum of completely irreducibles. It is enough 
to assume B is equal to a direct sum of completely irreducibles. Then, in addition to equality occurring 
in (3) when the right-hand side is 0, equality occurs precisely when the components of x are equal 
within each piece corresponding to a direct summand of B. 

For the third statement, it is enough to assume that B satisfies p (B) = 1. It then follows from the 
irreducibility of B that there is a positive diagonal matrix D such that DBD 1 is row stochastic (D~ l is 
obtained from the Perron-Frobenius eigenvector of B, which is positive). Since Bx = y if and only if 
(DBD 1 )Dx = Z)y, (3) holds for DBD' 1 if it holds for B, so that we may as well assume B is row 
stochastic. Then, if B is not doubly stochastic, some column sum is <1, so let = 3 6 (J < 1. Let 
X\ = 1, i =7^= /, and Xj = 1 -f~ c. Then (3) becomes ' 

1 + « < n (1 + M =l + $e + 0(e 2 ), 

i 

which is impossible if e > is small enough. This completes the proof of the theorem. 

It should be noted that essential portions of theorem 1 may also be demonstrated by a maximization 
argument. 

Example: The assumption of irreducibility in the third statement of theorem 1 cannot, in general, be 
relaxed. If B = ( i -i ) then B satisfies (3) for all x > and p (B) = 1, but B is not similar to a 
doubly stochastic matrix. 

An alternate form of (3) is 

Corollary 1: // B = (bjj) is an n by n doubly stochastic matrix, then for any n real numbers t\ 9 m * *,t n , 
satisfying U > — 1, i = 1, • • • ,n, we have 

n(l+t,) ^nd + stjbu). 

i = 1 i = 1 j = 1 

Our primary observation concerns row stochastic matrices with property (*) . 

Theorem 2 : If A is a row stochastic matrix with property ( * ) and D and E are the positive diagonal 
matrices guaranteed by ( 1 ) then 

(7) detDE>l. 

Furthermore, equality holds if and only if A is actually doubly stochastic. 

Proof: As B = (6y) runs through all n by n doubly stochastic matrices and F = diag {/ l9 - • -,/ n } 
runs through all positive diagonal matrices, then A = D _1 BF runs through all row stochastic matrices 

n n 

with property (*) where D = diag (3 b^fj,' ' ', 2 b nj fj}. Thus, since B = DAE, where E = F~ l , it 

suffices to show that det D > det F. If we denote (/i,- • -,/ w ) T by /, this is equivalent to saying that the 
product of the entries of Bf is greater than or equal to that of / for any positive vector /. This, of course 
follows from theorem 1. To analyze the case of equality in (7), it suffices to assume B is completely 
irreducible. In this event, it follows from theorem 1 and the fact that Bf has no components that 
equality in (7) implies that all entries of / are the same. Thus D = F and equality holds in (7) precisely 
when A is already doubly stochastic. 

Note: A related but rather different inequality when A is symmetric appears in [3, theorem 3]. Also a 
portion of the proof of that result could be used to prove part of the first statement in our theorem 1. 
It follows from theorem 2 that 

Corollary 2: If K is a row stochastic matrix with property (*) and B is related to A by (1), then 
|detA| < |detB|. 
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We denote the eigenvalues of A by ai," ■ *,a w , ordered so that |ai| < • • • < |a n |, and those of B by 
Pu m ' ',Pn, ordered so that |/?i| < • • • < |/? n [. Since a n = 1 = /? n , it follows from corollary 2 that 



Corollary 3: If A is a row stochastic matrix with property (*) <mrf B is related to A by (1), then 

n— 1 n— 1 

n |ai| < n |ft|. 



We conjecture that in case ^4 is row stochastic with property (*) and B is related to A by (1), 
then actually 

The result of theorem 2 may be extended to all matrices with property ( * ) in the following way. 

Theorem 3: // A is any nonnegative n by n matrix with property (*) a/i^ D a/irf E are positive 
diagonal matrices guaranteed by (1), then det DE > p(A)" n . Furthermore, equality holds only if 
DE = p(A) al L 

Proof: It is enough to assume A is completely irreducible (for, if not, it is equivalen to a direct 
sum of same) and then A is irreducible. In this event there is a positive vector x such that Ax = p{A)x 

and, therefore, ( X' l AX is row stochastic, where X = diag{#i, • • •,#»}. Application of theorem 2 to 

-^j7 If 1 AX yields det VE' > 1 where D'{-±- X 1 AX)E' = B. Setting D = ^tt D'Z" 1 and 
P (A) ^ - p (i4) p(^) 

£ = X£', gives fi = Z)^^ and det Z)^ > p{A)' n as was to be shown. The case of equality also follows 

from theorem 2. 

Remark 3: The reader may wish to note the relationship between the present work and the notion of the 
equilibrant, 

E(B) =ini P (FB) 

(where the inf is taken over all positive diagonal matrices of determinant 1), of a nonnegative matrix 
mentioned in [5] . 
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